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Abstract 

We analyze the Mott transition in multi-band Hubbard models with the inclusion of multiplet 
exchange splittings as it arises in infinite dimensions by using the generalized Gutzwiller wave- 
function introduced by Biinemann, Weber and Gebhard [Phys. Rev. B 57, 6896 (1998)]. We also 
present an extension of that variational wave-function to account for broken-symmetry solutions, 
which still allows an exact analytical treatment. Our analysis reveals some drawbacks of the 
variational wave-function, which, in our opinion, imply that Gutzwiller- type of wave-functions do 
not properly characterize quasi-particles close to a Mott transition. 
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I. INTRODUCTION 



After so many years since its proposaii, the Mott transition maintains intact its scientific 
interest, continually kept alive by the growing number of materials which, on the track of 
becoming Mott insulators, show rich and fascinating physical properties, among which high 
Tc superconductivity is likely the most spectacular example. 

Much of the theoretical effort has been till now devoted towards cuprate-inspired mainly 
single-band models, like the standard Hubbard model or its strongly correlated counterpart, 
the t — J model. Although unsolved questions still remain, a lot has been unveiled about 
the Mott transition in the single-band Hubbard model, especially after the development 
of the so-called Dynamical Mean Field Theory (DMFT)^. DMFT, which is exact in the 
limit of infinite coordination lattices, represents a quite reliable approach to investigate the 
on-site dynamical behavior across the Mott transition. Yet, in its original version, DMFT 
is not able to yield accurate results concerning inter-site correlations, which are treated in 
a mean-field like fashion, although several extensions has been proposed^. For that reason, 
other approaches has been often adopted to investigate the possible occurrence of d-wave 
superconductivity in proximity to the Mott insulating phase of single-band Hubbard and 
t — J models, most of which are more or less explicitly related to the Gutzwiller variational 
technique^. That amounts to study a variational wave-function consisting of a simple Slater 
determinant where doubly occupied sites are partially or completely projected out. The 
Gutzwiller wave-function can be rigorously handled only numerically in finite dimensions, 
although an approximate analytical scheme to evaluate average values was proposed by 
Gutzwiller himselS^, thereafter called Gutzwiller approximation formula (GAF). Later on, it 
was realized that the GAF is exact in the case of infinite coordination lattices^ and that it 
provides similar results to the slave-boson method within the saddle-point approximation^. 
Moreover, the comparison with the numerical treatment of the Gutzwiller wave-function has 
shown that the GAF is quite accurate also in finite dimensiona^iM. 

In spite of the interesting developments achieved by means of the Gutzwiller variational 
wave-function in single-band Hamiltonians, especially for what it concerns c/-wave super- 
conductivity in the t-J mode l^i^°i^^ , there have been not many attempts to apply the same 
technique to multi-band models, even though there are many interesting strongly correlated 
materials where orbital degrees of freedom play an important role. Indeed there exist several 
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analyzes based on multi-band extensions of the GAFi^ii^iii and of slave-boson techniquei^, 
but a direct comparison with more exact results to test the quality of the Gutzwiller wave- 
function is still missing, even though DMFT results for multi-band Hubbard models are by 
now availableiS^iiimo. 

In this paper, we review the multi-band extension of the Gutzwiller wave-function (GWF) 
proposed by Ref. Il^ which has the big advantage of being analytically treatable in the limit of 
infinite dimensions. Moreover, we propose a further extension to broken-symmetry phases, 
which maintains the same property of being analytically accessible in infinite dimensions 
without loosing any variational freedom. In this limit, we single out some peculiar properties 
of the GWF close to the Mott transition. The comparison with exact DMFT results reveals 
some drawbacks of the GWF which may lead to even qualitatively incorrect results. 

The plan of the paper is as follows. In Section |n] we present the generalized GWF 



introduced in Ref. 



121 for multi-band models and calculate the general expression of the 



variational energy in infinite dimensions. In Section III we study the infinite-dimension 
Mott transition both in the absence and in the presence of multiplet exchange splittings, 
which we discuss more in detail in Section IV. In Section V we present an extension of the 
GWF to symmetry broken phases which we use for a specific two-band model in Section VI. 
Finally, in Section VII we draw some conclusions. 



II. THE MODEL 



We consider a multi-band lattice model where, on each lattice site, N valence orbitals are 
available for being occupied by the conduction electrons. Since our purpose is to discuss some 
general features of the Gutzwiller wave-function, we will assume the most simple form of 
tight-binding energy with nearest neighbor hopping matrix elements diagonal in the orbital 
index, namely 

t ^ 

^0 = -7= mm 4,aaCj,aa + H.C., (l) 

V ^ (jj) a=l cr 

where z is the lattice coordination, Cj and c| are, respectively, the annihilation and 
creation operators for an electron at site i with orbital index a = 1, . . . ,N and spin a =t, i- 
The correlation among the electrons is introduced via two local interaction terms: an on-site 
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Hubbard repulsion 

Hu = ^Y.nl (2) 

i 

where = Ef=i E<t '^i aa'^i acF the electrou occupation number at site i, as well as an 
exchange splitting term 

2N 

Hj = ^^^Jv„\i,n,Vn){i,n,Vn\. (3) 
i ?i=o r„ 

Here |z,?T.,r„) denotes a multiplet of ra-electron states at site i. r„ has multiplicity g^,^ such 
that 

^ 2N 



H 9Tr. = 9n 

r„ 



V 



n 

being the binomial on the right hand side the total number Qn of available on-site n-electron 
states. Without loss of generality, we assume that Hj only splits multiplets at fixed n 
without affecting their center-of-gravity energy, implying I]r„5'r„<^r„ = 0. 

In the absence of interaction, the ground state I'I'o) of (0) with an average electron number 
per site no consists of degenerate bands each one at filling no/2N, namely 

l*0) = n n iaA,ai|0), (4) 

"=1 \k\<kp 

where 

V ~ N' 

\k\<kF 

If < rio < 2A^ is integer, we expect that, for sufficiently large U, the correlated ground 
state of H = Hq + Hu + Hj describes a Mott insulator. Analogously to what has been done 
for the single-band Hubbard model, we would like to have a system which allows to study 
a metal-to-insulator transition into an ideal paramagnetic Mott insulator, in contrast to a 
more conventional metal-to-band-insulator transition. For that reason, we further assume 
that the multiplet exchange splitting (jH)) is such that, within a perturbation expansion upon 
the uncorrelated ground state |$o); all self-energy diagrams are diagonal in spin and orbital 
indices and independent of them. This does not exclude that spin and/or orbital symmetry 
may be spontaneously broken especially close to the Mott transition, as it is known to occur 
for the single band Hubbard model on a bipartite lattice. Yet, in what follows, we will 
discard such a possibility and just discuss an idealized MIT, where neither the metal nor the 
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Mott insulator break any of the symmetries of the Hamiltonian. We postpone to Section 
the analysis of spontaneous symmetry breaking nearby the Mott transition. 

Having this in mind, we start analyzing the role of the interaction by a Gutzwiller vari- 
ational approach. Namely we search for the best variational wave-function of the form 

|^G)=nAG|'^>o), (5) 
i 

where the uncorrelated wave-function is the Slater determinant of Eq. (0]) and the on-site 
Gutzwiller projecto r^^i^"^ , 

Kg = ^^Kv„\i,n,Vn){i,n,Vn\, (6) 



tends to go along with the local interaction terms Hu and Hj in mo 



difying the relative 



12, there is a certain 



weights of on-site electronic configurations. As shown by Refs. Ig 
arbitrariness in the choice of the variational parameters A's, related to the fact that any 
transformation acting on |$o) and involving operators of which |$o) is an eigenstate amounts 
simply to a multiplicative factor. In our case, that arbitrariness allows to impose, without 
losing generality, the normalization condition 

{^o\PIg\^o), (7) 
as well as an additional constraint on the single-particle density matrix 

where the last equality derives from our choice of |$o), Eq- ©• 
A formal solution of (jZj) and (jH)) is obtained^ by writing 

P(o)(n,r„)' 

where P^^^n, F^) and P{n, r„) represent the occupation probabilities of the ra-electron mul- 
tiplet r„ in the uncorrelated, |$o), and correlated, I^E'g), wave-functions. Upon inserting (jH)) 
into ((Tj) and (jH} one obtains 

5:5:p(n,ro = i, (10) 

" r„ 

the correct normalization for P(n,r„), as well as 

^5:nP(n,r„) = no, (11) 



namely the condition that the average number of electrons per site coincides with the un- 
correlated value uq. In this representation, the correlated probability distribution is the 
variational quantity which has to be optimized. 

The Eqs. (j7)) and (jHI) imply that, within a perturbation expansion in the parameters 
(1 — A„r^)'s5 only more than two fermionic lines can exit from any vertex P^q at site i, 
see Ref. il2|. This property simplifies considerably the calculations in the limit of infinite 
coordination lattices, 2; ^ oo in Eq. (^, where one can show that the average variational 
energy per site^^ 

Eg = ^{^g\Ho + Hu + HjI'^g) 

= ZTo + ^Y.Y.i^-^^fnnJn) + Y.Y.Jr.P{n,V^). (12) 

To is the average value per site of the non-interacting tight-binding Hamiltonian 

To = lim h%\Ho\^o), 

which gets reduced by a factor Z through the Gutzwiller projection, and we have introduced 
a chemical potential term which makes the Hubbard interaction minimum at the average 
n = riQ. Our choice of Hj leads to the following expression 

It, — -L i jT, i-L ri— 1 

= E E J P{n,T^)P{n-l,T^-,). (13) 

n^i r„jt., 9n-i V no {2N - no) 

When Hj = 0, namely if nothing in the Hamiltonian splits the degeneracy among states 
at fixed n, there should exist a variational solution in which all configurations at fixed n are 
equally probable, a property owned by the Slater determinant |$o)- Iii that case, a multiplet 
with degeneracy gr„ occurs with probability 

P(n,rO = ^P(n), (14) 

9n 

being P{n) the probability of a site to be occupied by n electrons irrespective of the config- 
uration. Eq. ()14j) inserted into (jlHj) leads to the following simple expression of Z valid for 
Hj = 0: 

2^ n ' 



E / , J^^P{n)P{n - 1). (15) 

^ ^rio (2iV-no)V^?"-i 



n=l 
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III. MOTT TRANSITION WITHIN THE GUTZWILLER APPROACH 



The search for the optimal P{n, r)'s with arbitrary values of the interaction parameters 
is not straightforward. However it is still possible to obtain simple analytical results close to 
the Mott transition, which is signaled by a vanishing quasiparticle residue Z. Indeed, when 
Z <^ 1, one realizeai^i^i by inspection of Eq. (fT^ that 

P(n,r„)~^l"-"»l. (16) 

It is therefore justified to assume P{n, r„) = for all n's but n = tiq and n = uq ±1, which 
simplifies remarkably all calculations. 

In the simple case Hj = 0, which we denote hereafter as J = 0, this approximation leads 
through Eq. (fT3|) to 



no gno 



^riQ (2N - no) V 9^o-i 



P{no)P{no - 1) 



+ , """"^ ,^^P{no + l)P{no). (17) 

The distribution probabilities have to satisfy 

Pino - 1) = Pino + l) = d, 

and consequently 

Pino) = l-2d. 

By inserting those expressions into (fTTjl . one obtains 

Z = ^"^^ ^ f\/^o(2iV -no + l) + ^ino + l)(2iV - no)]' 

no(2A/-no) ^ J 

= jiN,no)dil-2d), (18) 

hence 

EG = 7iN,no)dil-2d)To + Ud. (19) 
The optimal d which minimizes fll9p is readily found: 

^ TojiN,no) + U ^ UciJ = 0)~U 

4To7(iV,no) " 4f/,( J = 0) ' ^ ^ 
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where 

f/e(J = 0) = -7(iV,no)To 



1 



2A^ - no L 



y'no(2iV -no + l) + v^K + l)(2Ar - no) 



no 
(21) 



is the value of the interaction at the Mott transitions^, when the optimal d = 0. The 
variational energy is therefore 



EG{J = 0) = 2Toj{N,no)d^ 



1 [u,{j = 0) - uy 



(22) 



8 Uc{J = 0) 

Let us now consider the case Hj ^ 0. We assume that the exchange splitting favors at 
any given n a particular multiplet of states, which we denote as F*. The Mott insulator 
described by the Gutzwiller wave-function is therefore characterized by P(no, r*^) = 1 and 
has energy Ems = Jr^^ < 0. In order to better understand how the Mott transition occurs 
when J 7^ 0, it is convenient to consider separately two extreme cases which do not require 
any numerical calculation. 

Let us start by considering an hypothetical metallic solution able to smoothly transform 
into the Mott insulator, in which therefore only the multiplets F* favored by Hj are occupied 
close to the MIT. Since only no± 1 and no are relevant, one obtains an expression of Z similar 
to ()18|) with the only difference that 



no(2A'— no 

I 1 2 

1 3^nQ+i 



v/no(2iV-no + l)./?^ 

V 9no-l 



9v , Qr* 

9no 



+ ^(no + l)(2iV-no)i 



9no+l 9no 



< l{N, no) 



The variational energy reads now 



Eg{J7^0) = ^,ff{N,no)d{l-2d)To + Ud 

+ ( Jr. , + Jr* - 2 Jr. ) d+Jr* , . 

V "0-1 "0+1 "0/ "0-1 

Provided we substitute 7 — > 7e// and 

U ^Ueff = U + f Jr. + Jr. - 2Jr. ) , 
^J-i V "0-1 "0+1 "0/ ' 

which is the actual Hubbard repulsion measured with respect to the energies of the lowest 
multiplets and not from the centers of gravity, the formal solution has the same expression as 
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before. In particular the critical interaction at which that metallic phase becomes unstable 
is now 



f/e( J ^ 0) = -7e//(iV, no) To - (Jr;^^_^ + Jr;^^^^ - 2 Jr.^ 

= J = 0) - (Jr. + Jr. - 2Jr. ) . (23) 

7(A^,no) ^ ' V ^no-i ^"0+1 ^"oy ^ > 

For very small J's, UJ^J 7^ 0) is shifted down with respect to VJ^J = 0) by terms of order 
|To|, which already suggests that the above solution is not the most energetically favorable, 
although it has the merit to merge smoothly into the insulator. 

Indeed, the best variational solution is actually different from the above one. If the J's 
are very small compared with To, we rather expect in the metallic phase that the probability 
distributions of the multiplets at any n are only slightly modified with respect to the J = 
case. If this were true, we could still search for a variational solution of the same form as for 
J = 0, which has an energy given by (j^UI) with d as in Eq. ()22|) . This solution, however, does 
not converge into the insulating one, which has an energy Ems = ^r*^^- The two energies 
indeed cross when 

1 [U,{J = 0)-Uf _ 

8 Uc{J = 0) ^"o' 

namely when U = U^, 



= U,{J = 0) - ^-8t/,(J = 0)Jr.^. (24) 

For small J, U^, is larger than Uc{J 7^ 0) given in Eq. (|^. which suggests not only that 
the J = metallic solution has lower energy but also that the Mott transition is first order. 
The explicit solution of the variational equations shows that, in the optimal metal at small 
J, the probability distributions are indeed only slightly modified with respect to the J = 



case by terms of order y | J/To| hence that the Mott transition becomes first order as soon 
finite J is introduced. 

In conclusion, if the multiplet exchange splitting term | J| is much smaller that the uncor- 
related bandwidth W, then the Gutzwiller variational approach leads to a first order phase 
transition from a metal, slightly modified with respect to the J = case, into a Mott insu- 
lator instead dominated by J. This transition is predicted to occur when the quasiparticle 
residue 

l^{N,no) (25) 



Tn 
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and arises because the metallic solution has to pay too much hopping-energy to modify the 
relative weights of the multiplets, a cost which overcomes the exchange splitting energy gain. 



A first order phase transition has indeed been found in Ref. |12 by explicit numerical 
minimization of the variational energy in a two-band model, and agrees with linearized- 
DMFT results obtained in Ref. 113 on a similar model. Moreover, in the same Ref. ilTj it is 
shown that the first order character reinforces with increasing exchange splitting strength J 
from J = 0, while it weakens with further increasing J above some intermediate value. This 
also agrees with our above results. Indeed, when J gets so large to make Uc{J 7^ 0), Eq. 
(EHjl . greater than f/*, Eq. (j21I), we expect the Mott transition to turn again into a second 
order one. That anyway requires a substantial \J\ ~ |ro|. 

IV. DRAWBACKS OF THE GUTZWILLER WAVE-FUNCTION 

In the previous Section we have identified the qualitative behavior across the Mott tran- 
sition for multi-band Hubbard models within the Gutzwiller variational approach. In par- 
ticular our analysis has been performed in the limit of infinite coordination lattices, where 
the Gutzwiller wave-function is able to describe a Mott transition. In the same limit, exact 
results can be obtained by DMFT— , which allows a direct comparison hence a test on the 
quality of the Gutzwiller wave-function. 

As we previously said, upon approaching the Mott transition with a Gutzwiller wave- 
function for multi-band models, the on-site charge configurations with a number of electrons 
n different from the average one uq get suppressed approximately like Zl""*^"', where Z is 
the quasiparticle residue. It is usually assumed that the GWF in the metallic phase of a 
single band model provides a faithful description of what Landau quasiparticles do, while it 
gives a very poor characterization of the underlying insulating component. If this property 
holds in multi-band models too, then the quasiparticle gas should be characterized by a 
P{n) ~ ^l"^"o|^ This behavior is not displayed by exact DMFT results, which rather 



indicate a quasiparticle probability distribution P(n) ~ Z (see Ref. Il6| ). The origin of this 
disagreement can be easily traced back. 

The main reason why the Gutzwiller wave-function gives a poor description of the Mott 
insulator in the single band Hubbard model at half-filling is that it can not account for virtual 
excitations into unfavorable charge configurations. Those virtual processes are responsible 
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of the finite average value of the double occupancy (n^ni) in the insulating phase even in 
infinite dimensions. Indeed {n-^ni) as function of U displays a discontinuity in the slope 
across the Mott transition in infinite dimensions^. The singular part, which vanishes at the 
MIT, is attributed to the quasiparticles and it is qualitatively reproduced by the behavior 
of the double occupancy as obtained through the Gutzwiller wave-function, supporting the 
common belief that this wave-function does correctly capture quasiparticle properties. 

In a multi-band Hubbard model the situation is different. If we just consider the occu- 
pation probability P{nQ ± 1) of {uq ± l)-charge configurations, we expect across the MIT a 
behavior similar to the double occupancy in the single-band Hubbard model, still compat- 
ible with the GWF. Let us instead consider the occupation probability P{nQ ±2). In the 
GWF close to the Mott transition, those configurations get suppressed like Z^. In reality, 
virtual processes from the more advantageous uq and (no ± 1) charge configurations imply 
first of all that P{nQ±2) is finite in the insulating phase too, and secondly that the singular 
(no ± 2) quasiparticle contribution still linearly vanishes across the MIT, unlike what it is 
found by the GWF. That disagreement is more profound than what it would seem to be. A 
quasiparticle probability distribution of the Gutzwiller type, namely P (n) ~ sug- 
gests that quasiparticles remain more strongly interacting than implied by the true behavior 
P{n) ~ Z, even after the Hubbard side-bands are well formed. This indicates that, unlike 
what happens in the single-band Hubbard model, the multi-band GWF is not fully adequate 
to capture quasiparticle properties. 

The second failure of the Gutzwiller variational approach regards the onset of the first 
order phase transition, which is predicted to occur when Z ~ \/\J]^ ioi \J\ <^ \Tq\. It also 
originates by the lack of virtual processes. If J = 0, the Gutzwiller wave-function leads 
to a Mott insulator with a finite entropy, related to the finite number of degenerate on-site 
electronic configurations with no electrons. This state has an infinite susceptibility to a term 
Hj which splits that degeneracy, with an energy gain linear in J. This result is obviously 
wrong. The super-exchange terms generated by virtual processes into unfavorable charge 
configurations lead to finite susceptibilities even in the Mott insulator. It implies that the 
actual energy gain is quadratic in J so that the Mott transition is either second order or 
weakly first order, in this case occurring when Z ~ \J\. 

Indeed, this aspect is not peculiar to a multi-band model but also occurs in a single 
band model in the presence of a magnetic field B which splits spin-up singly occupied sites 
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from spin-down ones. Also in that case the Gutzwiller approach would predict a first order 
transition when Z ~ while in reality, being the magnetic susceptibility finite^, the 

transition occurs at smaller Z, see Ref. likely when Z ~ 5. 

Yet those defects of the Gutzwiller wave-function might not affect qualitatively the physi- 
cal behavior in the most common situations where the multiplet exchange splitting term leads 
to the conventional Hund's rules, namely favors high spin and angular momentum configu- 
rations. These states are usually representable by a Slater determinant hence are accessible 
by a mean-field approach which can be improved by the Gutzwiller projector. However there 
may be less conventional but still interesting cases where the multiplet exchange splitting 
favors low degeneracy states, which are not Slater determinants hence unaccessible to mean 
field theories. Here the Gutzwiller wave-function might be inadequate to describe the Mott 
transition mainly because it is unable to access the interesting region where the metallic 
kinetic energy gain ~ -Z^ITqI competes with the exchange splitting ~ J. 

The inability of the GWF to describe phases which are not accessible by Hartree-Fock, 
like an ideal Mott insulator, is a well known fact in finite dimensional. However, it is com- 
mon opinion that the GWF is instead able to yield a correct description also for Hartree- 
Fock-unlike physical situations in infinite dimensions. Even more, it is believed that the 
infinite-dimensional picture provided by the GWF is qualitatively right even in finite di- 
mensions. The above discussion suggests that this common belief is partly wrong and that 
an improvement of the GWF towards including virtual inter-site processes is necessary^^. 

V. GUTZWILLER WAVE-FUNCTIONS FOR SYMMETRY BROKEN PHASES 

In spite of its appealing features, an ideal Mott insulator at zero temperature is unlikely to 
exist, especially if it has huge degeneracy. Commonly one expects a symmetry broken phase 
to occur at low temperature. For instance, in a single-band model at half-filling, the ideal 
Mott insulator has an infinite spin degeneracy which is likely reduced at low temperature by 
some magnetic ordering. Therefore, even though any mean-field type of approach (including 
more sophisticated ones based on Density Functional Theory) can only stabilize correlated 
insulators in broken-symmetry phases, namely can only describe band-insulators thus hiding 
the basic phenomena leading to a Mott insulator, yet they often provide a faithful description 
of the low temperature physics. 
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In this situation, the Gutzwiller variational approach should still be useful to improve 
Hartree-Fock approximation. That would amount to search for the best wave-function of 
the form 

\^g{A)) = Pg |$(Ao)) = n Ag I<^>(Ao)), (26) 

i 

with PiG still given by Eq. and where |$(Ao)) is a symmetry-broken uncorrelated 

trial wave-function with single-particle order parameter Aq. In general, after Gutzwiller 
projection, the correlated wave-function will have a different order parameter A. This implies 
that the average values of the single-particle density matrix over |\1/g'(A)) and |$(Ao)) do not 
coincide. This does not obviously represent a problem for a numerical treatment, whereas 
would seem to prevent the use of the method developed in Ref. 

the average values in infinite dimensions. In fact that method relies on the possibility of 
constructing a GWF with the same average value of the single particle density matrix as 
the uncorrelated wave-function. We could still impose that condition for the wave-function 
defined by (j^ . but that would reduce the variational freedom. 

In this Section we present a simple extension of the GWF to account for symmetry 
broken phases while leaving the property of being analytically treatable in infinite dimensions 
without any variational loss. We start by noticing that there always exists a non unitary 
operator U = HiUi such that its action 

t/|$(Ao)) = |$(A)), (27) 

leads to a trial wave-function |$(A)) of the same form as |$(Ao)) but with the same average 
value of the single-particle density matrix as the Gutzwiller projected |\I/g(A)), hence the 
same order parameter A. Therefore, 

|*g(A)) = Pg |$(Ao)) = Pg f> |$(A)) = Pg|$(A)), (28) 

which implies that, provided we substitute 

Pg^Pg = PgU, (29) 

we can still search without loss of generality for a variational wave-function 

|v1>g(A)) = Pg|$(A)), 
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where the average single-particle density matrix stays unchanged after Gutzwiller projection. 
The cost is that we must work with a Gutzwiller projector as given by (j^^ . which is in general 
neither diagonal in the multiplets which appears in Hj^ Eq. (jS)), nor hermitean. 

Actually we can identify two distinct situations which may occur, one of them being 
already included in the formalism developed by Ref. Q. If the exchange splitting J favors 
a degenerate atomic configuration, then we reasonably expect that in the true ground state 
of the lattice only one of the degenerate states will be occupied on a given site, eventually 
changing from site to site. There the order parameter corresponds locally to a conserved 
quantity of the atomic Hamiltonian, for instance the 2;-component of the on-site spin, which 
leads to a generalized Gutzwiller projector still hermitean and diagonal. Let us consider 
as a simple example a one band model. The local Gutzwiller projector is in general 

Pr,G = Ao K 0) (z, 0| + A2 K, 2) {i, 2| + Ai [ |z, T) (z, T I + i) i | ] , 

where |z,0(2)) denote the empty or doubly-occupied site i, while |2,cr) the singly occupied 
site with spin a. We assume that the uncorrelated wave-function is magnetically ordered. 
Then ipTjl may be constructed by local operators 

?7, = e^"»^S 

where Sf is the ^-component of the spin operator at site i. We find that 

Pi,G = Pi,G U^ = Ao \t, 0) {t, 0| + A2 \t, 2) {t, 2| + Ai e"' \t, T) {t,U + Ai e""' \t, i) 

Indeed the modified Gutzwiller projector is still diagonal and hermitean, although there 
appear different variational parameters for spin up and down components. This additional 
degree of freedom gets fixed once we require that the order parameters of the correlated and 
uncorrelated wave-functions coincide. The above type of GWF is actually a particular case 
of the generalized GWF's introduced in Ref. Q. 

We can however envisage a different situation where the uncorrelated wave- function |$o) 
has an order parameter which does not correspond locally to a conserved quantity. There 
a modified Gutzwiller projector Pq is unavoidably off-diagonal and non-hermitean. Let us 
discuss an over-simplified example. We assume that the exchange term leads to a local 
Gutzwiller projector of the form (we drop the site label) 

Pg = Xa\a){a\ + Xb\b){b\, 

14 



and moreover that the uncorrelated wave-function has an order parameter identified by 
non-zero matrix elements 

($o|a)(6|$o) = ($o|6)(a|$o). 
The transformation U can be now taken of the form 

U = a {\a){a\ + \h){h\)+ p {\a) {h\ + \h) {a\) . 

We readily find that 

Pg = PgU = a\a\a){a\ +a\b\b){h\ +/3AJa)(6| +/3Afe|6)(a|, 

indeed containing off-diagonal terms and evidently non-hermitean. This is a novel situation 
which we are going to discuss more in detail in a particular example. 

We conclude this Section by pointing out that the non-hermitean character plays a crucial 
role only if the ordering involves just the quasiparticles, while it is essentially irrelevant when 
both the quasiparticles and the Mott-Hubbard side-bands contribute to the order parameter. 
In the following section we analyse a two-band model where both cases may appear. 



VI. A TWO-BAND MODEL STUDY 

We consider a two-band Hubbard model described by the Hamiltonian and (j2)) where 
the orbital index a = 1,2. Besides the local spin-density operators 

^ I 2 t 

^ a=l aji 

where a^-, Oy and dz are the Pauli matrices, we introduce orbital pseudo-spin operators 

= 9 E E 4aa^a6C, (30) 

where the Pauli matrices f's act on the orbital indices. The hopping and Hubbard terms, 
and (0), have a very large SU(4) symmetry. Having in mind common physical realizations, 
like e.g. rf-orbitals of tg symmetry, we assume that SU(4) is lowered down to a spin SU(2) 
times an orbital 0(2) by the exchange term, which can therefore be written as 

Hj = Y.JsS,-S,^ JtT, .fi-3{Js + Jt) {Tff . (31) 
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The above Hj just splits the on-site configurations with two-electrons. There are six of 
those states on a given site i, a spin-triplet orbital-singlet \i,n — 2; S — 1, Sz;T — {)) , which 
we denote hereafter as \i,2,t), and a spin-singlet orbital-triplet, for which we use the short 

notations |i, n = 2; S* = 0; T = 1, = 0) = |i, 2, 0) and |i, n = 2; 5* = 0; T = 1, = ±1) = 
|i, 2, ±). In this subspace, Hj has the form 

Hj = Y^2Js \i, 2, t){i, 2, t| + 2 Jt \i, 2, 0)(^, 2, 0| - (3 J5 + Jt) K 2, ±)(^, 2, ±|. (32) 

i 

The standard Hund's rules correspond to — < Js < — 5 Jt/6 < 0, when the spin-triplet 
orbital-singlet configuration has the lowest energy, followed by the spin-singlet orbital dou- 
blet with Tz = ±1. In this case the ideal Mott insulator at half-filling, no = 2, represents 
localized spin-1 moments which should order at low enough temperature to freeze out the 
spin-entropy. On general grounds one expects that the magnetic ordering in the insulator 
should contaminate the nearby metallic phase so that, as U increases from weak coupling, 
first a transition from a paramagnetic into a magnetic metal should occur, followed by a 
Mott transition into a magnetically ordered insulator. Even a mean field approach is in 
principle able to reproduce the above scenario. In this situation, as we discussed before, 
the Gutzwiller-projected Hartree-Fock wave-function does improve the mean-field solution, 
providing a better physical description. In Appendix I we analyse in detail the case in which 
a bipartite lattice stabilizes an antiferromagnetic ordering. 

Less conventional is the situation Jt < —\Js\, where the non-degenerate spin-singlet with 
Tz = 0, namely 

\i, 2; 5 = 0; T = 1, 7; = 0) = {4,1^21 + 42t4u) 1°)' (33) 

is the lowest energy configuration. That would be for instance the case of two Hubbard 
models (two-chains, two-planes, etc.), coupled by an antiferromagnetic exchange. Here the 
large U Mott insulator with Uq = 2 describes a collection of on-site singlets, a local version 
of a valence-bond (VB) insulator. Since it is not degenerate and fully gaped, we expect the 
VB insulator to be stable at large U against any spin and/or orbital order. Just to avoid 
unessential complications, we assume that the lattice is sufficiently frustrated to prevent any 
spin/orbital ordering at any U. This situation is far less trivial than the previous one. In 
fact, being the Mott insulator not describable by a single Slater determinant, it is inherently 
unreachable by any mean-field approach, which necessarily leads to some kind of ordered 
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state. According to our previous discussion, we expect in this case that also the GWF is 
unable to provide a faithful description of the Mott transition. 

As we showed in Section 11111 the GWF without any symmetry breaking would undergo 
a first order metal-insulator transition when the quasiparticle residue 




namely when U = U^ being 

= -6To - A^GToJt. 

We also argued that this result is wrong since a metallic phase is able to enter the regime 
in which Z ~ \Jt/To\. Let us now check whether there exists a better broken-symmetry 
metallic solution. Indeed, even if lattice frustration prevents spin/orbital order, there is 
still a broken-symmetry GWF which might in principle compete with the above metallic 
solution. 

A. Superconducting Gutzwiller wave-function 

When U = 0, the multiplet exchange term (jH^ favors a BCS Hartree-Fock wave-function 
with the s-wave order parameter 

($Bcs(Ao)l4i^42j$Bcs(Ao)) = ($Bcs(Ao)l42t4iJ$BC5(Ao)) = Ao < i. (34) 

When Aq —>■ 1/2, the doubly occupied sites in the spin-triplet, \2,t), or in the doublet of 
spin-singlets, |2, ±), configurations get suppressed by a factor (1 — 2Ao)^ with respect to the 
Tz = spin-singlet, |2,0). Similarly, the probability of singly, |1), or triply-occupied, |3), 
sites vanish like (1 — 2Ao). This suggests that by Gutzwiller projecting out sites with zero 
and four electrons, |0) and |4), respectively, through the variational wave-function 

I^g(A)) =Pg|<^>bc5(Ao)), (35) 

one might indeed smoothly connect to the VB insulating state, with (1 — 2Ao) playing the 
role of Z. However, even though the uncorrelated BCS wave-function has a large order 
parameter Aq ~ 1/2, the correlated {"^g) should have a much smaller one, A ~ ZAq, since 
only quasiparticles are involved in superconductivity. In such a case we are therefore obliged 
to implement the non-hermitean Gutzwiller projector in order to get analytical results for 
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large coordination lattices. In other words, we shall work with a Gutzwiller wave-function of 
the same form as (jH^j) and impose that I^&g) and \^bcs) have the same order parameter A 
through a non-hermitean Pg, see Eq. (j29p . Since this is a novel situation in the Gutzwiller 
variational approach, we prefer to describe it in detail. 

In order to simplify the analysis at half-filling, we assume that the Hamiltonian has 
particle- hole symmetry and search for solutions which do not break it. This guarantees that 
there are still two conditions we can impose without losing variational freedom: an overall 
normalization, ((Tj), and a particle-hole symmetry constraint. 

To accomplish this job, it is convenient to work not in the original electron basis but in 
the natural basis where the on-site single-particle density matrix is diagonal. This is done 
by the following unitary transformation, which is valid at half-filling Hq = 2, 



( c 

«,1(2)T 
V 42(l)i 




^ '^i,l(2)T 



In the natural basis the only non-zero on-site average is 

(<l>BC5(A)|<,,a,,,,|<l>BC5(A)) = ^ - A. 

In order to distinguish the local configurations in the natural from the original basis we will 
denote the former ones as |n, r^j). The most general local Gutzwiller projector is in this case 

Pi,G = ^'^KTn \i-,n,Yn)^^%,n,Yf\ + Ao4 |i,0)(i,4| + A40 |2,4)(i,0|. (36) 

The last two terms are the only possible off-diagonal elements at half-filling when particle- 
hole symmetry holds. The normalization condition and the conservation of the single-particle 
density-matrix lead to the following parametrization of the A's for 7^ 0,4 



where Fj^) and P^^\n, F^) are the correlated and uncorrelated occupation probabilities 
in the natural basis. For n = 0, 4 we have instead 

A^P(°)(0) + A^4P(°)(4) = P(0), (37) 
AloP(°)(0) + AlpW(4) = P(4), (38) 
AoA4oP(°H0) + Ao4A4P(°)(4) = ^(0; 4) = A(4; 0), (39) 
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where we introduce the transition amphtudes 

The occupation probabihties in the natural basis are related to those in the original one 
through 



P(0) 


= 1 (P(0) + P(2, 0) + A{0; 4)) + v^A(0; 2, 0), 


P(l) 


= P(l) + A(l;3), 


P(2,0) 


= P(0)-A(0;4), 


P(2,±) 


= P{2,±), 


P(2,t) 


= P{2,t), 


P(3) 


= P(l)-A(l;3), 


P(4) 


= i (P(0) + P(2, 0) + A{0; 4)) - V2A{0; 2, 0), 


A(0;4) 


= |(P(0)-P(2,0)+A(0;4)), 



(40) 



where we have used the fact that, by particle-hole symmetry, P(l) = P(3) and P(0) = P(4). 
The order parameter A is given by 

2A = 2\/2.4(02+) + A(13), 

and the following inequalities should be verified 

2A(0;2,0)2 < P(2,0)[P(0) + A(0;4)], 
A{0;A) < P(0), 
A(l;3) < P(l). 

The big advantage in working with natural orbitals is that the Z reduction factor has the 
same expression as in (fT!?|l . namely 

T. ^«I^(l-A)""(l + A)"A,,,,A,_,r._,. (41) 

Once we know how to relate the parameters A's and Z to the variational occupation probabil- 
ities defining the correlated and uncorrelated wave-functions, we can solve the most general 
variational problem by minimizing the energy functional 

Eg{A) = Z{A) To(A) + U [4P(0) + P(l)] 
+2Js P(2, t) + 2Jt P(2, 0) - {3Js + Jt) P(2, ±), (42) 
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where 

{^Bcs{A)\Ho\^Bcs{A)) =n{A). 

For the sake of clarity, we present here an analysis based on the following parametrization 
of the A's in Eq. (I36|) . although it has less variational freedom: 



<(0) 



Ao - I (Ao + A20) Y -^^-^ 



A4 - ^ (Ao + A20) ]J -^r^y 
A - - \ ) 



'<(0) 



A40 - i (-^0 - A20) Y pm(o)' 

Ai =A,,/H, (43) 



A3 — All 



pr(i)' 
V pr(3)' 

A- - A 

^24 - A2a/-(o) 



A - A 

As before Aq is the order parameter of the uncorrelated BCS wave-function while A is the 
true order parameter after Gutzwiller projection, see Eq. (j^ . P^^{n,Tn) and P^\n,Tn) 
are the distribution probabilities in the natural basis for the BCS wave-functions with order 
parameter Aq and A, respectively. They are explicitly written in the Appendix, Eq. ()B4|) . 

The normalization condition as well as the conservation of the single particle density 
matrix imply that 

A„r„ - (44 j 

PAo (^, r„) 

where P^^{n,Tn) = ($AoNj r„)(z; n, r„|$Ao) is the occupation probability for configura- 
tions in the original electronic basis within the uncorrelated BCS wave-function with large 
order parameter Aq, see Eq. ()B5|1 . while P{n,Tn) is the same quantity for the correlated 
wave-function. Eq. ()44|) is the most natural generalization of (jH)) to a broken-symmetry 
phase, which is the reason why we have chosen the above parametrization. The true order 
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parameter A is defined through 



2 A = P^\0) + -Pf (1) - -Pr(3) + (4) 




+ Ao(l-4A2 




P(l) 



(45) 



which indeed is of order Z when 1/2 - Aq ~ P(2, 0) ~ 1, P(l) ~ Z and P(0) ~ Z^. The 
explicit evaluation of the Z reduction factor is presented in the Appendix, see Eq. ()B2|) . 

Let us now compare the variational energy as given by Eq. P2|l with Z of Eq. (jB6jl . 
valid for 1/2 — Aq = 5^1, to the energy of a non superconducting paramagnetic solution, 
Eq. (021) with A = 0, Z being given by Eq. ()B3j) . We find that the Gutzwiller projected 
BCS wave-function has always higher energy by terms roughly of order Z\Tq\. 

Therefore, even though the Gutzwiller projector is quite efficient to transform the huge 
Hartree-Fock energy cost, namely 



($Bcs(Ao)|i^|$Bcs(Ao)) - {^Bcsm^l'^Bcsm = To(Ao) - To(0) + 2UAl + UtAI 



into a much smaller one of order Z\To\ ~ ZU close to the Mott transition, yet it is not 
able to make superconductivity favorable. Namely, the best variational metallic solution 
remains the one described in Section IIIIl with all the drawbacks discussed in Section |TVl In 
conclusion, as we anticipated, the Gutzwiller variational approach does not properly describe 
a mean-field-unlike Mott transition. 

In reality we may expect a superconducting phase just before the VB Mott insulator. 
Recently it has been studied a model which shares many common features with the present 
one, namely a three-band Hubbard model with inverted Hund's rules^^, mimicking a strongly 
dynamical Jahn- Teller effect. For an average number of electrons per site uq = 2, the inverted 
Hund's rules favor, like in our example, a non degenerate singlet on-site configuration. By 
a DMFT calculation, a superconducting instability was discovered just before the singlet 
Mott insulator. However, that instability was found to appear when the quasiparticle residue 
Z ~ I J| (see Ref. Ilfil l. As we discussed at length previously, the simplest metallic Gutzwiller 
wave-function which we have so far considered is unable to reach Z ~ J, since it becomes 
disadvantageous with respect to the insulating one already at Z ~ Therefore we can not 
exclude that superconductivity may occur even in the two-band model we have considered. 



To{0) + - + Jt, 
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VII. CONCLUSIONS 



In this paper we have analysed some pecuhar features of the Mott transition displayed 
by a multi-band Gutzwiller variational wave-function (GWF) in infinite dimensions. The 

n 

analysis has been carried out by using the generalized GWF introduced in Ref. which 
allows a simple analytical treatment in infinite dimensions. Moreover, we have extended 
that wave-function to account for broken-symmetry phases. 

It is usually assumed that the GWF in infinite dimensions gives a faithful description 
of the quasiparticle behavior around the Mott metal-insulator transition. We have shown 
that while this belief is partly true for single-band models, it is incorrect for multi-band 
models. In particular, we have identified at least two major failures of the GWF across 
the Mott transition. The first concerns the occupation probability P{n) of on-site charge 
configurations with n electrons different from the integer average one no, which is believed 
to represent just the quasiparticle occupation probability normalized to the quasiparticle 
residue Z. The GWF in infinite dimensions predicts P(n) ~ ^l"~"o| close to the Mott 
transition, while both physical arguments as well as Dynamical Mean Field Theory results 
suggest a P{n ^ uq) ~ Z, even in infinite dimensions. This apparently innocuous disagree- 
ment is instead profound. In fact, the GWF results imply that the quasiparticles remain 
much more strongly interacting than what the correct P{n) ~ Z behavior suggests. 

Another drawback concerns the Mott transition in the presence of a weak multiplet 
exchange splitting term J. Within the GWF, the Mott transition turns into a first order 
one and occurs when the quasiparticle residue Z ~ yJ\J\/W, W being the bare bandwidth, 
much before the quasiparticle gas has had the time to react against J. This happens because 
the susceptibility to an infinitesimal exchange splitting J diverges at the Mott transition for a 
GWF. In reality, that susceptibility is finite so that the Mott transition is either second order 
or weakly first order, in that case occurring when Z ~ | J|/iy. The main consequence is that 
the interesting region where the metallic hopping-energy gain ~ ZW competes against the 
exchange J is not even accessible by a Gutzwiller wave-function. Both the above mentioned 
shortcomings have the same origin: the inability of the GWF to account for virtual processes 
into unfavorable charge configurations. 

We have then argued, on the basis of a two-band model study, that the GWF is still a good 
variational wave-function in all cases which can be qualitatively described by a mean-field 
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theory, but it fails otherwise, as for instance in the case we have exphcitly analysed where 
the Mott insulator is a local version of a valence bond insulator. There an improvement of 
the GWF is necessary. 
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APPENDIX A: ANTIFERROMAGNETIC GUTZWILLER WAVE-FUNCTION 
FOR A TWO-BAND MODEL 

The best Hartree-Fock wave-function |$o(^)) on a bipartite lattice when the multiplet 
exchange term favors the spin-triplet configuration is the ground state of an Hamiltonian 

Hhf = Ho + I^Y1 Y1 (-^T (^»,«T " ^im) > 

i a=l,2 

which describes an antiferromagnetic insulator with order parameter 

($o|rii,iT - ^imI^o) = (^o|rii,2T - ^i,2i|^o) = 2(-l)^m. 
We search for the optimal Gutzwiller wave-function 

|*G) = nAG|*o(m)), 

i 

where for a given sublattice and making use of particle-hole symmetry, 

Pi,G = Ao [\t,0){t,0\ + \tA){iA\] 

+Ai+ [ \i, 1; = 1/2) (z, 1; = 1/2| + |z, 3; = 1/2) (z, 3; = 1/2| ] 
+Ai_ [ \i, 1; = -l/2){i, 1; = -1/2| + \i, 3; = -1/2) {i, 3; = -1/2| ] 
+A2± \i, 2, ±)(i, 2, ±1 + A20 \i, 2, 0) {t, 2, 0| 
+A2t+ \t,2;S=l,S, = l)(^,2;5 = 1,S, = 1| 

+X2to \h 2; 5 = 1, 5, = 0)(i, 2-S^l,S,^ 0| 
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while for the other sublattice + with — interchange. In this particular case the Gutzwiller 
projector remains hermitean since the non-unitary transformation Ui is diagonal in the above 
multiplet basis. 

We define the correlated probability distributions for a given sublattice 

P(l, +) = P(l, = 1/2) = P(3, = 1/2), 

P(l, -) = P(l, = -1/2) = P(3, = -1/2), 

P(2,t+) = P(2,5=l,5, = 1), 

P(2,t-) = P(2,5, = -l), 

P(2,t0) = P(2,5, = 0), 

and analogously for the uncorrelated ones P^^\. . .). For the other sublattice 5*2 <-> —Sz- 
The two conditions ((Tj) and (jH)) imply that 



X2 _ 


P(0) 


P(0) 


P(o)(0) Q 


— 1 


A?+ = 






P(o)(l,+) 


2(| + ^)'G-^) 


A?- = 






P(o)(l,-) 


2(i-^)'(i + "^) 


\2 _ 


P(2,±) 


P(2,±) 


P(o)(2,±) 


2 (i - m^) 


X2 - 
^2,0 - 


P(2,0) 


P(2,0) 


P(o)(2,0) 




X2 - 
^2,t0 ~ 


P(2,t0) 


P(2,t0) 


P(o)(2,tO) 




X2 - 


P{2,t+) 
PW(2,t+) 


(i+m)" 


X^ - 


P(2,t-) 
P(o)(2,t-) 


P(2,t-) 



The correlated occupation probabilities satisfy the normalization condition 
2P(0) +2P(1, +) + 2P(1, -) +P(2, 0) +P(2, ±) +P(2, tO) + P(2, t+) + P(2, t-) = 1, (Al) 
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as well as the conservation of the order parameter 

P(l, +) + P(2, t+) - P(l, -) - P(2, t-) = 2m. (A2) 

By all the above defined quantities, the variational energy is found to be 

Ecim) = Z(m)To(m) + t/[4P(0) + P(l,+) + P(l,-)] + 2JTP(2,0) 

+2Js [P(2, t+) + P(2, to) + P(2, t-)] - {3Js + Jt) P(2, ±), (A3) 

where 

To(m) = (<l>o(m)|i^o|$oH), 

and 



Z(m) 



P(h+) P(h 



+ 2(Vni'+) + Vni'-)) IVn2,±) 



'P(2,0) , P(2,t0) 



^ / P(2,t+)P(1,+) ^ / P(2,t-)P(1, 



(A4) 

For any finite [/ the optimal solution has always m 7^ due to the nesting property. For very 
large U we expect m — ^ 1/2. In this limit we can neglect all P's but P(l, +) and P(2,t+) 
hence, from Eqs. (jAll) and (jA2|l . 



P(2,t+) = 4m - 1, P(l,+) = l-2m, 

which implies that 

, , 4m — 1 

Z(m) ~ 2 . 

^ ^ l + 2m 

In the same limit the uncorrelated hopping energy has the expression 



To(m) ~ -2^M2(2-4m), 

where 



M2 = I dep{e) e^ 



is the second moment of the uncorrelated density of states per spin and orbital, p(e). There- 
fore the variational energy as function of the order parameter m for U ^ \Tq\ is 

4m — 1 



EG{m) ~ -A-——JM2 (2 - 4m) + f/ (1 - 2m) + 2Js (4m - 1) , 
1 + 2m * 
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and it is optimized by 



leading to 



Mo 



m 



Eg ~ 2J; 



2 (U - Us) 
2M2 



2 ' 



U - AJ: 



s 



(A5) 



(A6) 



APPENDIX B: EVALUATION OF THE Z-FACTOR FOR THE TWO-BAND 
MODEL 



The explicit expression of the Z reduction factor in the two-band model of Section IVT 
allowing for a superconducting order parameter is, through Eq. (PT|) . 



- - A 



--Aj+-A,A,o^2+^ 



- - A 
2 



+ |a5A„ (i + a) (i - a) + (i + a) (i - a^ ' 
+ ^A3A,o(i + A) (i-A)%A3A,(l-A;' 

If we parametrize the A's according to Eq. PHj) and make use of we find 



1 



VI - 4A2 4 



P(0)P(1) 



^^;?(o)pr(i 



(0), 



>(0) 



+ 



<^(4)Pi:(3) 



(0), 



)(0), 



>(0), 



+ 



(0), 



,(0) 



<H2,o)pr;(3) 



(0), 



,(0), 



)(0)/ 



+ -v/P(2,0)P(l) 



pSi^)PZ!a 



{0)l 



^pi°)(2,o)prj(i) ■ ^pi:(2,o)p]rj(i 



>(0) 



+ 



<^(4)Pi:^(3) 



(0)/ 



)(0)/ 



){0), 



+ 



V2 



^v/P(2,t)P(l) + ^^P(2,±)P(l) 



+ 



(Bl) 



(B2) 



When A = Ao = 0, the above expression reduces to the Z-factor for a paramagnetic non 
superconducting solution, namely 



= 2JP(0)P(1) + V3A/^(l)P(2,t) 



+ V2JP{1)P{2,±) + JP{1)P{2,0) 



(B3) 
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The uncorrelated probabilities distributions in the natural basis with order parameter Aq 
(analogous expressions hold for A) are 

, 4 



<^(0) = (l + Ao) ^l-45 + 65^ 

= 4(i + Ao)'(i-Ao):.45-125^ 



Pi:')(2;0) = (l + Ao)'(i-Ao)'-5^ 
Pi°)(2; ±) = 2 (i + Ao)' (I - Ao)' ^ 25^ (B4) 
Pi°J(2;t) = 3(i + Ao)'(i-Ao)'^3<5^ 



<^(3) = 4 (i + Ao) (i - Ao) ^0, 
^^(4) = (l-Ao)'-O, 

where the last expressions on the left hand side correspond to the limit Aq = 1/2 — 5 with 
5^1. The uncorrelated occupation probabilities in the original electronic basis are readily 
obtained by the latter upon inverting Eq. (^UJ: 

<H0) = Pi::)(4) = (i + A^)'^i-5 + 25^ 

= Pi:^(3)=4(3^-A^)^25-65^ 

Pi°)(2, 0) = 3^ + |Ag + A^ ^ i - 25 + 36^ (B5) 

P^)(2,±) = 2(l-Ag)'^25^ 

In the limit of small 6, by inserting ()B4|) and ()B5|) into ()B2|) one finds at leading order 
(recalling that A ~ Z ^ 1) 



y/Z v^y^P(0)P(l)(l+5) + v/p(l)P(2,0) 



+ 



2/P(l)P(2,t)VP(l)P(2,±) 



1+5). (B6) 
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